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Abstract 

Iterated Bessel processes IV {t), t > 0, 7 > and their counterparts on 
hyperbolic spaces, i.e. hyperbolic Brownian motions B^^{t), t > are ex- 
amined and their probability laws derived. The higher-order partial differ- 
ential equations governing the distributions of lR{i) = \W {2R'' {t)), t > 
and jR{t) = i7?^(|27?^(t)|^), t > are obtained and discussed. Processes 
of the form R''{Tt), t > 0, B''P{Tt), t > where T* = inf{s : B{s) = t} 
are examined and numerous probability laws derived, including the Stu- 
dent law, the arcsin laws (also their asymmetric versions), the Lam- 
perti distribution of the ratio of independent positively skewed stable 
random variables and others. For the process R''{T^), t > Q { where 
= inf{s : B^{s) = t} and B'' is a Brownian motion with drift /x) 
the explicit probability law and the governing equation are obtained. For 
the hyperbolic Brownian motions on the Poincare half-spaces H2, we 
study B^^{Tt), t > and the corresponding governing equation. Iterated 
processes are useful in modelling motions of particles on fractures ideal- 
ized as Bessel processes (in Euclidean spaces) or as hyperbolic Brownian 
motions (in non-Euclidean spaces). 

Keywords: Bessel process, Modified Bessel functions. Hyperbolic Brownian 
motions, Fox functions. Student distribution, Subordinators 
AMS: Primary 60J65, 60J60, 26A33. 

1 Introduction 

The analysis of the composition of different types of stochastic processes has 
recently received a certain attention with the publication of a series of papers 
(see for example [3], [2], [1], [7]). The prototype of these composed processes is 
the iterated Brownian motion whose investigation was started in the middle of 
the '90s. Beside the distributional properties of the composed processes much 
work was done in order to derive the equations governing their probability laws. 
It was found that these processes are related both to fractional equations and to 
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higher-order equations as is the case of iterated Brownian motion. The core of 
this paper considers the Bessel process t > started at point x G [0, oo) 

and with parameter 7 > at different random times. We firstly study the 
process = i?7(^2 ^ > where Rj, i?2 are independent Bessel processes 
with the same parameter 7. This is equivalent to studying the Bessel process 
at a random time which is represented by an independent Bessel process. 
Iterated processes have proved to be suitable for describing the motion of gas 
particles in cracks (or fractures). For the iterated Brownian motion this is 
considered in DcBlassie's paper [B] but a similar interpretation can be given to 
iterated processes obtained by composing Bessel processes (this is the case here) 
or fractional Brownian motions (see [7]). The law of Iii{t), t > is expressed in 
terms of Fox functions and possesses a Mellin transform equal to 



E{iR{t)r-' = {2H)^ ' ' 4 ' , 3?M > 1-7, t > (1.1) 

-.2) 



with 7 > and for /i = 7n + 1 produces the m-th order moments of Ifiit)^ t > 

-p / 7+m \ p / 27+m \ 

E{In{t)r = {2h)^^^^l^^y^, m = l,2,.... (1.2) 

We are able to present the p.d.c. satisfied by the distribution of lR{t), t > 
which differs for 7 > 1 and 7 < 1 because in the latter case an impulse delta 
function appears as in the iterated Brownian motion. The equation we obtained 
reduces to the fourth-order heat-equation 

dq _ 1 , 1 dH{x) 



dt 2^dx^ 2V2^ dx^ 

for 7 = 1. A related process considered in Section [5] is 

jR{t) = RjilR'm^), t>Q (1.4) 

where Rj, Rj arc independent Bessel processes starting at the origin. The 
probability density of (jl.4p can be expressed in closed form as 

q{r,t) ^ Pr{Mt) £ dr}/dr ^ ^.^^J^J^^^^ K^ (^) . r,t>Q (1.5) 

where Ki^ is the modified Bessel function of order zero (see [HI formula 3.478]). 
The equation corresponding to (|1.5p has the form 

dq l/a^a d^q h^lfdq - if \ , ^ „ 

and includes the equations governing the process Bi[\B2{t)\^) , i > 0, for 7 = 1 
(and coincides with 3.16 of [7] for = 1/2). Interesting results can be obtained 
by considering the Bessel process R'^{t), t > stopped at the first-passage time 
Tt, t > of an independent Brownian motion. Processes stopped at different 
types of random times can be viewed as processes with a new clock which 
is regulated by an independent Brownian motion B. The r.v.'s Tt = inf{s : 
B{s) = t} tells the time at which the Bessel process must be examined. This 
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means that the clock considered below is timed by an independent Brownian 
motion. Therefore W{Tt), t > Q represents a motion where accelerations and 
decelerations of time occur continuously. We show that the distribution function 
oi R?[Tt),t>0 reads 



Pr{W{Tt) G dr}/dr = 



1 



We show also that 
Pr ' ^ 



e dr\ /dr = ^ ^ 



1 



r,t>0. 



rA>0 



(1.7) 



(1.8) 



which coincides for 7 = ?i, f = 1 
Related distributions are 



Pr 



1 



£ dr> /dr 



n with the folded Student distribution. 

(1.9) 



^ r(i) 



Furthermore, we obtain also that 



t 

i+t2 



Pr 



t2 + |i?7(rop " J 



£ dr > /dr = 



_t2 



< w < 1. 



1 - 



t 



(1.10) 



with < r < t and 7 > 0. For 7 = 1 we have the arcsin law. Besscl processes 
Ry{Tl^), t > 0, 7 > stopped at first passage times Tt = inf{s : B^{s) = 
t} where _B'^ is Brownian motion with drift /i arc examined in Section 3. In 
particular we prove that 



q^{r,t) ^Pr{R^Tt) E dr} /dr 



4t e*^rT- 



2ir(i) 



i2 



i2 



(1.11) 



with r > 0, t > 0, 7 > 0, At > 0. 

The last section is devoted to compositions involving the hyperbolic Brown- 
ian motion, that is a diffusion on the Poincare upper half-space i?" = {xi, . . . , x„ : 
Xn > 0} with particular attention to the planar case and the three-dimensional 
Poincare half-space H'^. In the space we study the hyperbolic distance from 
the origin of a hyperbolic Brownian motion stopped at the first-passage time Tt, 
t > of the standard Brownian motion whose probability law can be explicitly 
written as 



sinh?] 



•y/cosh Lp — cosh rj 



7K0 



\^/¥Tp^]d^ (1.12) 



with 77 > 0, t > 0. In the distribution of J^{t), t > reads 
2\/2 7/ 1 sinh rj 



PJ3 i) 



TT (772 +t2)" 



(1.13) 
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The evaluation of the integrals leading to (|1.12p and (|1.13p necessitates the 
following formula 



exp {-l3xP - jx-P} dx = ^ (^^^ (2v^) (1.14) 



where p,j,l3,u > and is the modified Bessel function (see [HI formula 
3.478]). The equations governing p.l2p and (|1.13p are respectively 

d^p.js _ d^pj^ d ( 1 



dt^^ d,^ -^d-,[M^H^ .>0,i>0. (1.16) 

The hyperbolic distance of a hyperbolic Brownian motion plays in the non- 
Euclidean spaces if", n = 2,3,... the same role of Bessel processes in the 
Euclidean spaces. The structure of the probability law of B2^{t), t > is 
rather complicated (see formulae (j4.9p and (|4.1ip below) and therefore we have 
restricted ourselves only to compositions involving first-passage times. Much 
more flexibility is allowed by three-dimensional hyperbolic Brownian motion 
B^P{t), t > 0. Millson formula (see [10]), in principle, permits us to examine 
compositions of higher-dimensional hyperbolic Brownian motions stopped at 
random times. 



2 Composition of Bessel Processes with differ- 
ent types of processes 

We first present some information about the Bessel process of order 7 > and 
starting from a; > 0. The Bessel process K].{t), < > is a diffusion with law 



p(t,r;0,x) 



r ( r\2' 
1 \x) 



exp 



2t 



governed by the infinitesimal generator 



2 1 dx'^ X dx 



The Bessel function Iu{z) is defined as 



, a-,r > 0, t > (2.1) 



(2.2) 



(z/2)'- 



-2k 



k=0 



fc!r(fc + l^+l)' 



(2.3) 



For 7 = n , n S N, the Bessel process represents the Euclidean distance of a 
n-dimensional Brownian motion {Bi{t), . . . , Bn{t)), t > from the origin. For 



\ 



(2.4) 
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the explicit law (|2.ip reads 



xr 



p{t,r;0,x)^-[^-j- exp<j __i>/„_i(_), x,r > 0, t > (2.5) 

and simplifies for = as 



n — l 

p(t,r;0,0)=2-^-^^^ =r"-ifc(r,i), r > 0, i > (2.6) 



where the function k ~ k{r,t), r E M+. i > is the heat kernel satisfying the 
p.d.e. 



dt 2 ^ dr \ Or _ 

where A is that in (|2.2p with 7 = n. We point out that the probability density 
(|2.6p satisfies the p.d.e. 

9p 1 fc^^p , ^. d /'p\^ 1 fc^^p 7 — 1 9p 7— 1 1 

' > ^ ' -p^ (2.8) 



9i 2 [ dr\rJj 2 [ Qr^ r dr j-2 ^ 

for 7 = ?i. The differential operator figuring in (|2.8p is usually referred to as the 
adjoint of p.2p . For the operators introduced above the following interesting 
fact turns out to be very useful 

^=A*p and ^^Ak. (2.9) 

Our interest here is to study the composition of the Bessel process iR2{t), 
t > Q outlined above with different processes or also with an independent Bessel 
process 2Ro{t), t > 0. 

2.1 The Iterated Bessel process 

We consider here the iterated Bessel process 

iRit) ^ ,R^oi2RUt)), t>0 (2.10) 

where iRq, 2^0 independent Bessel processes of dimension 7 starting at 
a: = and possessing density 

2 

7— 1 — — 

p(r, t) = 2 ^ '\ , r > 0, t > 0. (2.11) 

^ ' (2t)*r(f) - ^ ' 

The distribution of (PTTU)) reads 



q{r,t) = 2' I —^-^—^-j-^ds, r>0,t>0 (2.12) 

,2; 



(2.)ir(2 (2t)irm 



and has Mellin transform equal to 

M{q{;t)}{r,) = (23t)^£li±iz^£ii±i_il, ^{,,} > 1-7. (2.13) 
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The distribution (|2.12p can be expressed in terms of Fox functions which are 
defined in the following manner: 



fj™^"/™"! TTm,n 

p,q \ ) p,q 

where 6* £ M and 



(ai,ai)i=i,..,p 



1 
27ri 



>-\-ioo 



M'^:^iv)x-^dr, (2.14) 



(2.15) 



By direct inspection of (j2.13p we see that in our case n = 0, m = 2, p = q = 2, 
bi ^ 7/2-1/2, 62 = 7/2-1/4, /3i = 1/2, /?2 = 1/4, ai = a2 = 7/2, ai = ^2 = 0. 
By considering the following property of the Mcllin transform 



(2.16) 



we conclude that (j2.12p can be represented in terms of Fox functions as 



(23t)l/4 2,2 



(23^)1/4 



(i,0); (i,0) 

_ 1 iV ('2 _ 1 1 

V2 2'2/' V2 4'4 



In view of the property of Fox functions (see [15] ) 



(6j +C^j,/3j)j = l,..,q 



c e 



for c = 1 , ()2.17p can be written as 



r(23t)i/4 2,2 



(23^)1/4 



(i,0); (i,0) 

V 2 ' 2-'' V 2 ' 4'' 



(2.17) 



(2.18) 



(2.19) 



Moreover, the Mellin transform p.l3[) generalizes the moments of (|2.10p and 
after some easy calculation, we have 



E{in{t)r = 2^ ^ ^ ^' m 



r2(|) 



m = 1,2, 



(2.20) 



Formula p.20|) shows that the mean distance of the iterated Bessel process 
increases as and the dimension 7 of the space where the iterated Bessel 
processes develop enters only in the multiplying coefficient in p.20p . In the 
next theorem we derive the governing equation of the law (|2.12p of the process 

Theorem 2.1. The density junction (|2.12p satisfies, forO < 7 < 1 the equation 



dq 1 f 1) ^ M f^' I 7-1 g I (7-l)(3-27) 



dt 23 V ar2 
1 



dr 



2{2t)iT {^) \dr^ ^^arrj'^^'^^- 



(2.21) 



For 7 > 1 the governing equation becomes 

dq _l_ f 8P_ 
at ^ 23 [d^ 



^ = 1 _ (^ _ 1)^1^ + A + (7-l)(3-27) ^ 



dr \ 



q. (2.22) 
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Proof. It is convenient to write the distribution (|2.12p as 

poo 

q{r,t) = / p{r,s)p{s,t)ds 
Jo 

where 



(2.23) 



p(r,t) = 2- 



'(2t)ir(i) r>0, i>0, 7>0. 

We start our proof by evaluating the time-derivative of p.23|) as 



dq 
dt 



f°° d 
- I p{r, s)—p{s,t)ds = [in view of (|2.8p ] 
Jo 

f°° 1 f 92 Q i-) 



7-1 



ds 



p{s,t) 



1 

~2 

We note that 



d_ 

ds' 



d 



P{r,s) <{ -Q^Pis,t) - 



(7-1) 



p{s, t) > ds 



= (5(r) hm = — - — ^ 



(2.24) 



and thus for s — > 0^ and s ^ 00 vanishes for all 7 > 0. A further integration 
by parts yields that 



d d d 

^P('',s)^p(s,i)f^s =^P(?', s)p{s,t) 



We note that 



— p(7-, s)p{s,t) 



i=0 



d 

= r'^~^-g^Hr,s)p{s,t) 



s=0 



i=0 



00 Q2 

, ^P{r,s)p{s,t)ds. 
ds^ 



[by (23] 



for < 7 < 1, because 



''"I / 92 7 - 1 a \ 5{r) 



2 ydr"^ r dr J r^-^ {2t)iT {^) 



(2.25) 



ds^"-'"' 2 

and p(r, s) — >■ (5(r) for s — > 0^. Furthermore, 

2 

p{s,t) = const s~'~^e~^ 

and this explains why 7 < 1 implies the appearance of the delta function. By 
collecting all pieces together we have that 



dq _r''-i f , 7 - 1 9 ^^ 5{r) 
'dt ~~ 



1 



9r2 r dr J r^~^2t)iT {^) 



(2.26) 
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By exploiting the following manipulations 

1 d r^^^ d 

^j,7-i / 52 d\l d p{r,s) 

2 \(?r-^ r dr J r dr rT~i 



equation (|2.26p takes the form 



dq ^ f & 



32 



dt 23 V9r2 ^ ; ^7-1 

7- 1 a \ 1 a 





r 








-1 


f 92 


r(i) 





22 \^ Qj.2 J, dr J r dr r'' ^ 

1 r-'-i / 92 7 - 1 9 \ (5(r) 



(2.27) 
(2.28) 



dr J r'T" 



In order to transform equation (|2.27p into the form of the statement of Theorem 
12.11 we note that 

^ - '7Z2-(7-l)7r:-U (2.29) 



^ 9r2 r dr J r^ ^ rT ^ \ dr^ dr r 

as a direct check shows. Furthermore, if we write 
92 d 1\ 

9;:2 -(T-i)a;-j'?(r-,i) = ^i'(^,i) 

and apply again (j2.29p we arrive at the following explicit expression 

92 dV q 1/92 , a 

-(7-l)Tr:- 9- 



c^T*2 7^ (^T* / 7.7 1 1 \ dr^ dr r 

By applying the same trick as above we have that 
52 7 — q 



Qj,2 J, dr J r dr r'^ ^ 
92 7 - 1 a \ 1 r (7-1) Id 



rT 1 \ dr"^ dr r J \ r^ r dr 

An analogous step must be applied to the singular term involving the Dirac 
delta function so that result (|2.2ip follows. □ 

Remark 2.1. For 7 = 1 equation (j2.2ip becomes 

dq _ 1 d*q 1 (P5{r) 
'dt~¥d?^^ 2V2^ dr^ 
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which is the fourth-order equation governing the iterated Brownian motion. 
This is because the process t > for 7 = 1 becomes the reflected iterated 

Brownian motion. It is weh-known that the law of the iterated Brownian motion 
satisfies the fractional equation 

diq 1 d'^q 

Remark 2.2. For 7 = |, equation (j2.22l) takes the form 

dt ~23 yOr"^ ^ 2drr) ydr"^ ^ 2r dr ) 

"23 ydr'^ ~ 22r2 c'r2 ^ 2r^dr)'^~ 2^ ^ 2^ dr dr J J 

This is the simplest equation involving iterated Bessel functions. 
If we consider the slightly modified process 

Jfl(t)== ii?;j(|2i?2(t)n, t>Q (2.30) 

we have the advantage that the law can be expressed explicitly in terms of 
modified Bessel functions. In fact, we have that the probability density g of 
jR{t), t > reads 



27.7-1 r°° I ^2 ^ 22r''~^ / 
g(r,t) = , — / -e-^-^ds = . —Ko[^\ (2.31) 

'^'^ (40^r2(2)io s (4t)^r2(2) "yvtj ^ 



where 



f°° 1 ( 1 
Kq{x) — J - cxp — I ds, x>0 

(see formula 3.478]). 

Theorem 2.2. The probability law of the process jR{t), t > is governed by 
the following third-order equation 

dq If d^q „,„ ,92(7 9(7 (7— 1)2] ^ ,„ „„^ 

Proof. We start by working on the probability density (j2.3ip which can be 
written as 

q{r,t)^ / p{r,s)l{s,t)ds, r>0,t>0 (2.33) 
Jo 

with 



p{r,s)^2-—^——, r>0, s>0 (2.34) 
(2s)2r(§) 



and 



l{s, t) = — — , s > 0, t > (2.35) 
(2i)ir(^) ^ ^ 
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We note that the density (|2.35p is a solution of the following equation 



dl „ dH , ^,dl 



(2.36) 



as a direct check shows (see also [TS])- In view of (j2.36p from (|2.35D we have 
that 



dq 
dt 



p{r, s) 



ds 



-2J^ s—p{r,s)—l{s,t)ds~{-f-2)J^ p(r,s)—lis,t)ds 

.00 Q 



= - 2 



°° d d 

s-^p{r,s)—l{s,t)ds + {-f-2) 



ds 



p{r, s)l{s, t)ds 



+2 



si d 



p{r, s) 



(7>0) 



7-1 9 
r dr 



22 yg^2 



I fc(r, s)l{s, t)ds 



We observe that 



92 


7 - 


1 9 


Qj.2 


r 


dr 


92 


7 - 


1 9 




r 


9r 


f 92 


+ ^ 


- 1 


I 9^2 




r 



s k{r, s)l{s, t)ds 



r 

sJ 



) fc(r, s)l{s, t)ds 



7— l9\/9 7— 1 



9r / \ dr 



and 



9 ^ 7-1 
9r r 



-'lir,t) = ^[l + r-]qir,t) 



dr 



Thus, we obtain 



dq 
dt 



^7-1 f g2 ^ ^ - 1 a 

Qj,2 



1 



r dr J r'^ 

2 [dr^-^^-'^d;--y^'^'^ 
1 / 9- 



9 



9 1 



by (1123] = -7T hzo - (7 - l)7r:r 1 + 9(^0 



7 92 , ^ 9 1\ , , 



2 V9r2 
9 1 
or r 



9r r 



9 



9?' 
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1 / 5^ 

2 [dr^ 



(7-1) d I (7-1) 
r (9r 



^^-(7-1^ 



g(r,i) 



Formula (|2.32p can be derived by simple calculation. 



□ 



2.2 The Bessel process at first-passage times 

Let Tt ~ inf{s : B{s) ~ t} where i? is a Brownian motion (possibly with drift 
yu) independent from the Bessel process ^{t), t > starting from zero. In this 
section we study the new process W{Tt), t > concentrating our attention 
on its law and sonic related distributions. Stopping the Bessel process W at 
the random time Tt can cause cither a slowing down (with respect to the nat- 
ural time) or a speed up of the time flow. The probability of slowing down is 
measured by the following integral 



Pr{Tt < t} 



te 2x 
/o VZTra?^ 

which decreases for all t since 

at V27r 
Furthermore, we observe that 

Pr{Tt <t] ' ^ 



dx 



t > 0. 



t\/2nx 



-.dx < 



e 2 



(2.37) 



(2.38) 



(2.39) 



and this confirms the asymptotic speed up of the time flow implied by the 
subordinator Tt, t > 0. We have now the explicit distribution of W{Tt), t > Q. 



Theorem 2.3. The distribution of R?{Tt), t > reads 



q(r,t) = P{W{Tt) e dr}/dr = 2 



r(i±i) 



tr'i- 



V^rli) (r2 + t2)' 



r,t>0. (2.40) 



Proof. 



q{r,t) =2 



r'^ 2s te 2s 



(2s)ir(^)V2^ 
2^^7-1 



zds 



2^0Fr,2 

2tr''-^ 

2^r^-ir(^) 1 

r(2) (r^+t^)"^ 







1 /r^+t^ 



e ™w 2 '^dw 



r,t>0. 



□ 
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Remark 2.3. We check that the distribution (|2.40p integrates to unity, 



^ r(^) J, (^2 + ^2)^ 



dr 



[r = tVv] 



1 r 



fr±l\ r°° 



^ r(i) 7o (l + y)^ 



1 + 2/ 



1 



+ 1 / w 



dw 



1 ~ w J (1 — U))2 



1 r 



«;2 1(1 — ?«) 2 irfiy — 1. 



With similar calculations wc can obtain the /i-moments for < /i < 1 



r''Pr{R'^{Tt) £ dr} - ^ ^ 



^ r(i) 



r(i)' 



\ — w J (1 — it;)2 
(l-u.)' 



V^r(i) 



We note also that for 7 = 1, the density law (|2.40p coincides with a folded 
Cauchy with scale parameter t and location parameter equal to zero. 



For the distribution function of R''{Tt), t > wc have the following result 

-dy (2.41) 



Pr{R''{Tt)>r}^ ^ ^ 2 



00 yy-l 



^ r(i) J, (^2+y2) 



2t r(^) 



r7-2 



7-2 



7-3 



—dy 



(7-l)(t2 + r2)^ (t2+y2) 

for 7 > 1. The recursive formula (|2.4ip yields some interesting special cases 

7 = 2, 

(f - arctan f ) , 7 = 3, (2.42) 
7 = 4. 



(42+^2)1/2 ■ 

Pr{R^{T,)>r}={ f ' 



t2+r2 
t 



((2+^2)1/2 ^ 2(t2+r2)3/2 

An interesting related distribution is presented in the next theorem. 
Theorem 2.4. The process 
1 



R^t) 



l + R^iTtY 
has distribution 

Pr{W{t) e dw}/dw 



t > 



(2.43) 
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r(i) ^ ^ 



t 

l + t2 



If 



l + t2 



< w < 1. 

(2.44) 



Proof. 



Pr{iV{t) e 

/ Pr{W{Tt) e dr} 
aw 1 1-'" 



(2.45) 



2< r(3ii) I n-w 



Hi r(^) U.2 



7-1 



1 — ii; 



2t r 



^7+1' 



-(1-u;) 



7-1 



2^ rf^) 



1 



\/rTt^; (i + t2)2 



w 



< w < 1. 



□ 



Remark 2.4. For 7 = 1 the distribution (|2.44l) offers the foUowing expression 



Pr{R^{t) e dw} =- 



t 



A 



(2.46) 



We are able to check that (|2.46p integrates to unity. Indeed 



A 



dw ~- 



TT Jq (w- 5)2+^2 TT J_B 1+y 



dy 



= — < arctan ■ 

TT 



1 - B 



A 



B 

arctan — 
A 



2 A 

— — arctan — —5 = 1 

TT A^ +B^ - B 

because A^ + B^ — B ~ {). In the general case we can verify that (|2.44p integrates 
to unity for all 7 > in the following manner 

/•l 7 /"OO 

Pr{IV{t) edw}= — Pr{R-^{Tt) e dr} 

.In dw 1 1--W 



w^ 
I — w 
w 



Pr \ R^{Tt) e d 



1 ~ w 



w 



poo 

/ Pr{R^Tt) e dy} = 1. 
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Remark 2.5. As a by-product of our calculations we show that 

l-T 



t 



1 +i2 



1 



-B 



1 i 
2' 2 



(2.47) 



We start from the relationship in the proof of Theorem 12.41 integrated in (0, 1) 



2i r(^) 1 (I 



r (^) 
" r(i) 



w 



7-1 



t2 



7-1 



t 



Now put 
obtain 



t 



i+t^) + l+t2 

^ = J/ in the first integral and w 



= ill the second and 



dy 



Thus the identity 

J-i (l + y2) — 
emerges. Therefore, 



4^ V 1 + 12 



-1(1 + ?/)^ 



77-1 



i (l + y2) 



+ y^ 



2r(2±i) 



Remark 2.6. Another result related to distribution (|2.40[) states that 



Pr 



t^ + \Ry{Tt) 



e dr }> /dr = - ^ 



1 - - 
t 



r\2- 



(2.48) 



(2.49) 



for < r < t, 7 > 0. It sufficies to evaluate the following integral 



Pr 



t^ + \Ry{Tt) 



2t 



—dy- 



r(i) drj^ 

For 7 = 1 from (|2.49p one obtains the law of sojourn time on (0, 00) of Brown- 
ian motion and the even-order pseudo processes, while for odd values of 7 the 
distribution of the sojourn time on the half-line for odd-order pseudo processes 
emerges (for 7 = 3 see [16], 7 = 2?! + 1, n > 2 see [TT]'). 

Remark 2.7. For 7 = n the process R"{Tt), t > can be represented as 



i?"(Tt; 



(2.50) 



where Bj{t), t > 0, j = 1, 2, . . . , n are independent Brownian motions and the 
r.v.'s Bj{Tt), t > possess Cauchy distribution. Therefore p.50p represents the 
Euclidean distance of an n-dimensional Cauchy random vector ( Ci (t ) , . . . , C2 (t ) ) , 
t > 0. 
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Remark 2.8. For us it is relevant to study the distribution of 

Sit) = — t > 0. 
After some calculation we find that 

P{S{t)edr}^dr^-^- — — ^, r,t>0. (2.51) 

We note that for t = , j = n the density (|2.5ip coincides with a folded 
i-distribution with n degrees of freedom and its density takes the form 

/(r;n)= ^ ^, r > 0. (2.52) 

For n = 1, the density (|2.52p coincides with a folded Cauchy and coincides with 
for 7 = 1 and at time t = 1. 

Theorem 2.5. The probability law (|2.40p o/ the process R?{Tt), t > is a 
solution to the following equation 

^ ' -(7-l)^-U, r,t>0, 7>0. (2.53) 



dt'^ \ dr^ dr r ^ 

Proof. It is easy to check that the density of the first-passage time 

t 

/(s,t) = -^, s,t>0 (2.54) 
satisfies the equation 

^ = 2J-, M>0. (2.55) 
In view of (j2.55p we have that 



--2p{r,s)fis,t) 



dr^ dr r 



d 

^P{r,s)f{s,t)ds 
OS 



□ 

Remark 2.9. For 7 = 1 the Bessel process coincides with the reflected Brown- 
ian motion so that R^{Tt), t > is a reflected Brownian motion stopped at the 
random time Tt and therefore becomes a folded Cauchy process. It is easy to 
prove that the Cauchy density 

q{r.t)= J r,t>0 (2.56) 

n[t'' + r^) 

solves the Laplace equation and this agrees with (|2.53p . 
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By inverting the role of the Bessel process and of the first-passage time we 
obtain a new process somehow related to W{Tt), t>0 which we denote by 

= inf{s : B{s) = W{t)} (2.57) 

where TfjT(j), i > is the first instant where a Brownian motion B first attains 
the level R^{t), i > (and K* is independent from B). The probability law of 
(|2.57p becomes 



2 



/ s'^e °2 [l + \)ds 
Jo 



A(x + t)* r(i) 

It can be checked that the distribution (j2.58[) integrates to unity. 

Remark 2.10. It can be easily checked that the following relationship holds 



r^.(t2) =■ R^Tt), t > 0. (2.59) 

From p.59p one can also infer that 

Tn-,it)=- {R-'{T^))\ t>0. (2.60) 
In particular, for 7=1 the result (|2.60p says that 

Tn^t)'^- {BiT^,)f '^- (ciVt))' , t>0. (2.61) 

Remark 2.11. Wc note that the probability density (j2.58p for 7 = 1, t = 1 
coincides with the ratio of two independent first-passage times through level 1 
of two independent Brownian motions. In other words we have that 

Pr \Tri(i) e dw] = Pr \Wi/2 G dw] = -^^^diy, u; > (2.62) 
L w J I / n w + 1 

where = '^1/9/^1/2 ^^^'^ -^1/2' ^1/2 ''^^'^ first-passage times of B^ and 

B^ through level 1 and are stable r.v.'s of order 1/2. 

Remark 2.12. The last statement is a special case of the following result. For 
stable positive, independent r.v.'s Tj, with Laplace transform 

Ee^^'^" = e^^" , A > 0, < < 1 (2.63) 
the ratio Tj/T^ where 

Pr{W^ e dw] /dw = ^ , u; > (2.64) 

(see e.g. [5], [H])- We provide a simple and self-contained proof of this result 
based on Mellin transforms. Let us consider two independent, positively skewed 
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stable r.v.'s Yi, I2 of degree v > with Laplace transform (|2.63p . The density 
function g of the ratio Y1/Y2 reads 



g{w) = / xf{x)f{xw)dx, w > (2.65) 
Jo 

where / is the density of Yi and I2 ■ The Mcllin transform of (|2.65p becomes 

/•oo poo f 1*00 ^ 

{Mg){'q) = / w""'^^ g{w)dw ^ / w"^^^ I / xf{x)f{xw)dx\dw 



^{Mg)iv)iMg){2-^^), < 3?M < 1. (2.66) 
In order to write / we resort to the Fourier transform (see e.g. [12]) and write 

f{^) = ^ cxp [^t/3x - |/3r cos ^ (1 - zsgn(/3) tan ^) } dp. (2.67) 

This is because 

Ee'^^ = [by (mSD] = exp{-(-i/3)''} exp {-|/3|''e-''^^sn(«| 

= cxp{-|/3rcos^ (l-isgn(/3)tan^)}. (2.68) 



We now evaluate 



1/31"" exp {-z^sgn(/3) - |/3re-^^^sn(0) | d/3 
exp {-z^ - /J-^e-^^ } d/3 

/•OO 

J /3-"exp|i^ -/3''e^'^|d;3 



Since 



/■oo /"OO 

I P~"exp[-iY-P''e~''^]dp + J /3-''exp|^^ -/3''e''^|d^ 







/3-''cxp{-/3''e-*'^}d/3 + e''¥ / /?"'' exp {-/J^e*"^ } d/3 

"'0 



1/ \ V ) V V ^ 

2 / 7r\ /1-?A 2 , , /I-77 
— cos TTTi — — r = — sm VKT]) r 



we have 



(7W/)(,j) = sin (TTjy) T {^—^ ) , < SRjr;} < 1. (2.69) 
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Analogously 



{Mf){2 -T]) = ^ — ^ sin (27r - ttt]) T 



77-1 



(2.70) 



r(2 - 77) 



sin (ttj^) r 



,7-i\ r(i-77) 



and thus the Mellin transform of g becomes 

sin^ (ttt]) r 



sin (Tr?^) r 



7] — 1 



77-1 



and in light of r(a;)r(l — x) = tt/ sinTra: we obtain 

1-77 



{Mg)irj)=—Tl^ ^ 



r 1- 



sin (7777) 



sin (7r77) 



v sm TT 



Let us take the density 

X sin nv 



h{x) 



n(x^ + 1 + 2x cos j/tt) 
and evaluate its Mellin transform 



sm 



, x>0, 0<i/<l 



dx 



x'^ + 1 + 2x cos T^TT 



(2.71) 



(2.72) 



(2.73) 



(2.74) 



We have that 

x'^'^h{x)dx 



smTTZ^ 



TT Jo {x + e-''-^''){x + e^^'') 







1 

1 

'27ri 

_ r(7; + l) 
27ri 

r(7; + l)r(-77) 



1 



1 



dx 







(x + e-*'^'') [x + e^^") 

poo 











dx 



27ri ^ 

r(77 + l)r(— 77) . sin7rj^7/ sin7ri'77 

sm7ri^77 = : = — : . 

TT sin— 777r siu 777r 



Furthermore 



TT^x^" + 1 + 2a;'' cos niy) 



x'^-'^h{x'')dx = - I y'^-'h{y)dy = 



dx, ^{t]} > 1 - 



1 sin TTT] — TT 1 sin 7r77 



sm - — TT I' sm -TT 



where the change of variable has been introduced in the first step. 
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Remark 2.13. A relationship between Mittag-LefHer functions and the distri- 
bution appearing in (j2.64p exists and reads 



e x'^ ^ sin niy 

/o 7r(a;2'^ + 1 + 2a:'^cos7ri/) 

with A>0, 0<i^<l. 



dx = Ei,,,{-Xf) (2.75) 



We now consider the composition Z(t) = C{Si,{t)), t > where C is a 
Cauchy process independent from the stable law Si, with Laplace transform 

^g-AS„(t) ^ g-tA-^ A > 0, < < 1, t > 0. (2.76) 

We remark that for u = 1/2, (|2.76p gives Ee~^^^/^'^*'> — e~*^ which represents 
the Laplace transform of the first-passage time of Brownian motion for level 
t/y/2, t > 0. The probability distribution of Z{t), t > reads 

f°° t 

Pr{Z{t) <E dz}/dz ^ — —p^{s,t)ds 
Jo + 5 j 

1 /•+00 i<ao 

_ / ^-^Px-m / p^{^s^t)dsdp (2.77) 



27r „ „ 

where p^{s,t), s > 0, t > is the law corresponding to (|2.76p . Therefore, in 
view of (j2.76p we have that 

1 f°° 

Pr{Z{t)edz}/dz^— e-"^''-*^'^^"dl3, xeE.,t>0. (2.78) 

3 Some generalized compositions 

We somehow generalize the previous results by considering the twice iterated 
Brownian first-passage time. By 

Ir^(t) = T^2 = inf{si : B\si) = inf{s2 : B^{s2) = t}}, t > (3.1) 

we mean a process which represents the first instant where a Brownian 
motion hits the level and represents the first instant where a Brownian 
motion B^ (independent of B^) hits level t. Clearly the probability density of 
the process (|3.ip reads 



Pr{lTit)edx}/dx^ -—=—=ds. (3.2) 

Jo V /TTX'^ V ZTTS"^ 

For the n-stagc iterated Brownian first-passage time we have 

/^(<)=Ti_i, t>0 (3.3) 

and the corresponding law becomes 

/"(t, x) =Pr{/^^ e dx}/dx (3.4) 
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, .... , — ^ ... , dsi . . . dSn-\ 

For the ri-timcs iterated Brownian first-passage time 

I^{t) = inf{si : B\si) - inf{s2 : ^^(sa) = . . . = inf{s„ : B(s„) = t}}} (3.5) 
we have the following theorem. 

Theorem 3.1. The distribution p.4p o/ the process T^{t), t > satisfies the 
following p.d.e. 

^rit,,) = 2'^-'^r{t,x) (3.6) 

and possesses Laplace transform 

e^^^/"(t, x)dx = exp 2^-^ } . (3.7) 

Proof. In view of (j2.55p . by successive integrations by parts we have that 



r\t,x)^ (3.8) 

dsi . . . dSn-l = 



dt^ 



V2^ ' ■ ' j2TTsl^, 



°° sie-^ S2e ^„_i 9^" ie 



Jo V27rx3 i9s2"_i J2ttsI_^ 



The Laplace transform of p.4p has a very nice structure and reads 

e-^^/"(i,x)dx = (3.9) 



oo 



e / . . . / . . . dsi . . . dsn-idx 

Iq V27ra;^ \2ttsI^^ 



^s^V2\S2e 2=1 te 2=->-i 

e , ^ ■ ■ ■ , asi . . . dsn-i = 

^2nsl_, 

e-*^27^ = exp{~iA^2i-^}, i > m > 0. 
If we take the Laplace transform of equation (|3.8p we get that 

^£(t,A)-22"-V^(i,A). (3.10) 
It is straightforward to realize that p.9p satisfies equation (|3.10p . □ 
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The 2" initial conditions pertaining to p.6p are given in the following form 



dx 



t=o 



(3.11) 



Remark 3.1. If wc consider the generalization of (|2.40p . that is 
the corresponding probability law satisfies the 2"-th order equation 



(9^ q 



d 



d 1 



ar^ or r ' 



(3.12) 



Remark 3.2. The following shows that between the iterated first-passage time 
It and the iterated Brownian motion there is a strict connection. Indeed the 
distribution (13.21) can be written as 



Pr{lT{t) e dx}/dx 



se 2x ie 2s t 
=ds = — 



V27ra3 V27rs3 x 



'27rx \'2ns 



zds 



-Pr{B'WB^{x)\) e dt}/dt, x>0, t>0. 



We now consider some subordinated processes involving the first-passage 
time of a Brownian motion with drift fi, say T/*, t > 0. To make this topic 
as self-contained as possible we present a derivation of the distribution of T/', 
t > 0. 

Lemma 3.1. The maximal distribution of a Brownian motion with drift fi reads 



Pr { max Bn 

Q<s<t ^ 



{s) > A ^Pr {t^ < t} , /3 > 0, t > 



(3.13) 



e 2 



Vt 



dw 



Vt 



dw. 



Proof. For /3 > a we have that 



Pr I max B^(s) > /3, B^{t) < a| =Pr {t^ < t, B,,{t) < a} 



(3.14) 



^-E Pr [P^it) < a\B^iT^)]] 



Furthermore, 



< a 
dw 



B. 



: exp ■ 



dw e2^('"-'3) 



: exp ■ 



(3.15) 



iw~p-f,{t-T^))^' 
2{t - 

{w-P + ^^jt-T^f 

2{t - n) 
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2 



2{t - T^) 



= [2f]-w^y]^ / , exp ^ 

Thus by inserting p.l5[) into p.l4p we get that 

Pr\ max B^(s) > /3, B^it) < a 

In the hght of the calculation above we can write that 



Pr I max Bn(s) > B 

' 0<s<t ^ 



Pr { max B^{s) > /3, B^{t) <p}+Pr {B^it) > /?} 

I 0<s<i J 



e 



e 2t 



g2p/3 / + / '^^dw = Pr {t^' < t} 



This concludes the proof. □ 
We list here some consequences of Lemma 13.11 : 

i) the density of reads 

qil3,t)=Pr{T^ edt}/dt = (3 ^ , /3>0, t>0 (3.16) 

ii) the density q{/3,t) satisfies equation 

iii) the following relations hold 







p-/3|Ai| ( Ap-2;3|Ml „ < 



r'P{r^^ e dt} = ^^-^ j if,_i/2(H/3), (3.19) 

/•oo 

/ e-^'Prm e dt} = e/^Me-^*\M+^ (3 20) 

Jo 

for /3 > 0, /i e M, ?7 > i A > 0. 
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We now present the distribution of W{T^), t > Q which generahzes the resuU 

Theorem 3.2. The distribution of ^{T^), t > Q has the following form 
q^{r,t) ^Pr{W{Tt^) G dr] /dr 



with r > 0, t > 0, /X e M. 
Proof. 



K 



"1+1 



•2+^2^ (3.21) 



1,ir,t)= —^—-Pr{Tteds} (3.22) 

Jo (2s)2r(^) 



2t e^^r-^- 



1 



2i^(^)^/2^ 



'e 2s 2 (Is 



(^) ' ^^'^ (ImIv/^^) . (3.23) 



4te*''rT-i / ^2 
~2ir(^) V2^ 

Result p.2ip emerges on applying formula (|1.14p . □ 



Remark 3.3. By applying the asymptotic formula for the Modified Bessel 
function 

Kv{x)^ for a;^0+ (3.24) 



(see p 136 [T^ or p. 929 [9|) we have that 

4 tr-^-^ 2^r(2±l) 

'^''^'■'^^"2ir(2)V2^(r2+i2)4i (r2+t2)4^ 



tr^-^ /7 + l\ 



r(i)0F(r2 + t2)^ V 2 

The equation governing the distribution (r, t) is given in the next theorem. 

Theorem 3.3. The distribution p.22p solves the following p.d.e. 
9 92 \ / Q2 g i\ 

Proof. Wc apply the derivatives — to distribution 

qf,{r,t)^ / p(r,s)p^(.s,t)ds, r > 0, t > 0, ^ > 0. (3.26) 



We readily have that (in light of (j3.17p ) 



dp „ / 92 Q l\ 

^'^l Ts^'^'^P'^^'^'^ds = TO] = - - (7- 9.(r,i). 

□ 

For /i = in p.25p one retrieves equation (|2.53p . 
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4 Compositions of hyperbolic Brownian motions 
on the Poincare half-space 

We consider the classical model of hyperbolic space represented by the Poincare 
half-space H2 = {x, y : y > 0} with the metric 



y2 



(4.1) 



The hyperbolic Brownian motion on H2 is the diffusion process with generator 
H defined as 



2 \ dx'^ dy^ 

and its transition function is the solution to the Cauchy problem 

^ = n2U, xeR, y>0 
u{x,y,0) = S{y - l)S{x). 



(4.2) 



(4.3) 



It is convenient to study the hyperbolic Brownian motion in terms of hyperbolic 
coordinates (77, a) where rj is the hyperbolic distance of (x, y) from the origin 
(0, 1) of H2 ■ In explicit terms (??, a) and {x,y) are related by 



cosh 1] 



x^ +1/2 + 1 

2^ ■ 



Furthermore, a is connected with {x, y) by 

.2 ,j^2_i 



tanct ~ 



2x 



We note that the formulas transforming (x, y) into (77, a) are 




sinh rj cos a 
cos 77— sinh r] sin a ' 

cosh sinh r; sin a ' 



77 > 



i<«<f 



(4.4) 



(4.5) 



(4.6) 



Some details on these formulas can be found in [17], [1]. The Cauchy problem 
(|4.3|) can be converted into hyperbolic coordinates as follows 



du 
'dt 



sinh rj dij \ ^ drj J sinh^ 77 da'^ 



77 > 0, t > (4.7) 



subject to the initial condition 

1/(77, a, 0) = (5(7/) for all ae[0, 27r). 

If we concentrate on the distribution of the hyperbolic distance of the Brownian 
motion particle from the origin we disregard the dependence in (j4.7p from a 
and study 



2 sinh 7] drj \ ' dii J 

u{ri, 0) — 5(77). 



du 
dt 



(4.8) 
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It is well-known that the solution to (|4.8p has the following form 

g(77, t) = ^ / " =d(p, 77 > 0, <> 0. (4.9) 

V7ri3 J „ vcosh V3 - cosh 77 



If we pass from problem (|4.8p to 

at ^ sinh77a>, y'^^^^^'ldvj " (4 10) 

(by means of the change of variable t' = t/2) we obtain a somewhat different 
distribution which reads 

t' 

— ^ /.OO — r 

g(7?, i') = / , r ■ 77 > 0, i > (4.11) 

^7r(2t)3 7,, vcosh - cosh 77 

In the first paper [S] (and also in the subsequent literature) the factor 1/2 does 
not appear and the heat kernel is (|4.11D (up to some constants). Of course the 
probability density p2(r],t) is given by 

P2{r/,t) ^ ainh 7] q{r],t), 77 > 0, t > 0. (4.12) 

A detailed derivation of (j4.9p and (|4.1ip is given in [T3]. 
We note that the probability distribution (|4.12p 

sinh77e~8 f°° we~^t , , ^ 

P2iv,t) = 7=^/ , i , --d^, 77>0, i>0 (4.13) 

VTTt^ Jn V cosh cp - cosh 77 



solves the adjoint equation 
dp2 1 ( d'^p2 d 



(coth?7P2)^, 77>0, t>0 (4.14) 



dt 2 [ 9772 a?7 

For the distribution (|4.13|) further characterizations are possible. Indeed, we 
can rewrite P2('7j t) as follows 

Mv,t)==V2e-iE[l,r...,^=^^], 77 > 0, i > (4.15) 

where the mean is taken with respect to the distribution of Tt, t > 0, which is 
the first-passage time of standard Brownian motion. Moreover, 

where R^{t), t > is the 2-dimensional Bessel process described above and the 
mean- value is taken with respect to the distribution of Bessel process . 
We give now an alternative form of p2(j],t) in terms of the Euclidean distance. 
Indeed, the distribution of B^P{t), i > in H2 can be written as 

e~8 d [°° ipe~^ I 

^2(^7, ^) = — 2—^— / a/ cosh 95 — cosh 7; d(p (4-17) 
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e » d 

= - 2- 



TTt drj 



{t}>ri] V cosh K^{t) — cosh 77 1 , T] > 0, t > 0. 



If we take an Euchdean right triangle with one cathetus of length cosh rj and 
hypotenuse y/coship, then cosh ip — cosh 77 represents the length of the sec- 
ond cathetus. Therefore the integrals above represent the length of the second 
cathetus weighted by means of the probability distribution of the Bessel pro- 
cess in the plane. Thus, the formula (j4.17p hightlights the relation between the 
distribution of the hyperbolic distance in H2 and the corresponding Euclidean 
distance in R^. We can recognize the additional factor of (|4.17p as a gamma 
distribution with parameters 1/2, 1/8. 

We now examine the hyperbolic Brownian motion B^P{t), t > Q stopped at 
the first-passage time Tt, t > of an independent standard Brownian motion B 
defined as Tt = inf{s : B{s) = t}. In other words we study the process 

J2{t)^B^PiTt), t>0 (4.18) 

with distribution 

PJ2iv,t)=J^ ^^^''''^''TI^'^'*' V>0,t>0 (4-19) 



-t sinhr/ / / r^dsdp 

Jo Jtj vcoshp - coshi] s^y'2TT 

t sinh rj f°° p dp 1 t- f ^ 



V2^n Jri Vcosh(p - cosh 77 {p^ + t'^) 



K2 Qv^^^^ 



where we have used formula (|1.14l) . The functions are related by the follow- 
ing recursive formulas 

K,+i{x) = K,^i{x) + 2-K,{x) (4.20) 

X 

and thus 

K2{x)^Ko{x) + -Ki{x). (4.21) 

X 

The function is the so-called Modified Bessel function of the second order. 
In analogy with the representation (|4.17p we can give the following expression 
for the distribution of hyperbolic Brownian motion stopped at Tt, t > 

pj,(77,t) = -^l^i? |c(i)Vcosha - cosh 77X2 (^^{C{t)f + t^^ \c(t)>n] 

with 77 > 0, t > 0, where C{t), t > is a Cauchy process. For the distribution 
(|4.19p we can state the following result. 

Theorem 4.1. The process (j4.18D has distribution (j4.19p . say pj^ = pj2(riit) 
which solves the following Cauchy problem 

?PJ2 = hr^-7r7n- 9(7/, o) = 5(77), r,,t>o. (4.22) 



dt-^ ^ \ drf drj tanh 77 
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Proof. The distribution of (|4.18p can be written as 



where ^(.s, t) is the distribution of Tf , t > 0. In view of (|4.14l) we have therefore 
that 

Q2 /-oo Q 

-^PJ2=-/ P2{v,s)2—g{s,t)ds 

f°° d d f 1 

=2 / g{s, t)—p2{v, s)ds = —^pj, - — —7-pj, 



ds ' drf' ^ drj ytanh?; 



□ 



The hyperboUc Brownian motion on ~ {x, y, z; z > 0} is the diffusion 
with generator 

2 f a^ ^2 g2 ^ Q 
\ dx^ ^ dy"^ ^ dz'^ \ dz ^ ^ 

The distribution of the hyperbohc distance of Brownian motion in possesses 
the form 

Smb. Ti G ^ 

P3(ry,t) =^3(77,*) sinh^ 7y = r ??e~~, 77 > 0, t > 0. (4.24) 

where k^{rj^t) is the heat kerneL The kernel 

^3(7?, t) - / "^.^^ , 77 > 0, ^ > (4.25) 
2V7rt3 smh77 

is the solution to 



dt sinh"^ T) di] \^ ' dri 



(4.26) 



By means of the transformation t = t' /2 equation (j4.26p is converted into 

du 1 1 d f . ^2 d\ , , 

— = - 5— — sinh77— M 4.27 

dt 2sinh2 77 977V dv J 

and formula (|4.24p takes the form 

, , „sinh77e~2 ,,2 sinh?? „ „ 

P3{v, t) =2 ' 7?e-^ = -e-'^^siv, t), 77 > 0, t > 0. 4.28 

V27rt3 V 

The law qsirj^t) is the distribution of the 3-dimensional Bessel process or the 
Euclidean distance on M^. Alternative form for the distribution ^3(77,7;) involve 
the first-passage time of a standard Brownian motion 

t smh 71 

P3{r],t) =2 sinh 776" 2 Pr{T^ e dt}/dt = 2—^Pr{T^ e dt}/dt (4.29) 
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where T,, = mf{s : B{s) = 77}, = mf{s : B^{s) = t]} with denoting 
a Brownian motion with drift equal to 1 and B a Brownian motion without 
drift. It is wch-known that the probabihty law qs^rj^t) appearing in (j4.28p can 
be expressed in terms of distribution of one dimensional standard Brownian 
motion B{t), t > since 

q3{'q,t)df] = 2Pr{maxB{s) E dr]} - Pr{B{t) e di]}, -q > 0, t > 0. (4.30) 

s 

From ()4.28p we can obtain that 

^ f sinhi?3(i) 1 t „ 

because on integrating w.r.t. 7/ both members of (|4.28p we obtain 

1= p3{il,t)dTj ^ e-i ^^q^(rj,t)dr] (4.32) 
Jo Jo V 

where q-i{r],t) is the distribution of R^{t), t > which is 3-dimensional Bessel 
process. Moreover, distribution (|4.24p solves the p.d.e. 

dp3 d'^ps „ 9 / 1 



dt drf df] \tanh77?7 ^'^J ^ 
which is the adjoint of the operator appearing in (j4.26p . 

The process Jslt), t > obtained by composing the 3-dimensional hyperbolic 
Brownian motion B^^{t), t > with Tt ~ inf{s : B{s) — t} with B independent 
from Bt^P, has a probability law equal to 

ie 2= _ _ f e is te 2s 



PJsiv,*) ^ P3iv,s)—^==ds = rj sinhi] e — — ^^==ds 
Jo vzTTS'^ Jo 2y/TTS2 vzns-^ 

For the governing equation of (j4.34p we present the following result. 

Theorem 4.2. The distribution of J3{t), t>Q, say pj^ = PJ^iVi^)) solves 

d'^ / d I \ 

-^PJ3=h5Z2- 21^173:3 b^a, rj,t>Q (4.35) 



dt'^ ^ \ drf drj tanh rj 
subject to the initial condition 

PjAv,0) = S{v). (4.36) 
Proof. Distribution (|4.34[) can be written as 



PJsiv,*) = / P3{v)9is,t)ds 
Jo 

where g{s,t) is the distribution of Tt, t > 0. We have 

"^^'"'3=-/ P3{v,s)2—g{s,t)ds 

f°° d / d"^ d I 



□ 
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